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Model considerations are always very important for the concept clarification in 
theoretical and experimental physics. In particular, an exact solvable model 
may provide much more detailed informations which are useful for a complete 
understanding of the nature of the phenomena. Thus, in the study of magnetic 
' phase transitions of inhomogeneous systems, different ways of modeling the sys- 

, tems can be found in the literature. But to implement the inhomogeneity of 

■ real systems into a model usually make the model itself hardly to have exact so- 
lution. However, inhomogeneous systems are frequently encountered in nature, 

^ ' and they also play an important role in material science. For example, to de- 

^ velope new materials, the inhomogeneous systems formed by doping with other 

materials in a bulk system are often grown artificially in laboratory. Thence, 
it is important to construct some exact solvable models so that we can deepen 
our insights to the nature of such systems. The hierarchical spin models may 
be viewed as one of the efforts in this direction. 
^ ' About one and half decades ago, spin models defined on hierarchical lat- 

, tices received much attention in the literature[l — 9]. Hierarchical lattices are 

defined as the infinite limit of iterative steps of replacing each fundamental 
bond or cell with a given type of decoration unit, and thence they are self- 

■ similar fractal lattices. Since hierarchical lattices are highly inhomogeneous in 
, the local connectivity and lack of translational invariance, spin models defined 

on hierarchical lattices may provide a very good frame in understanding the 
magnetic properties of inhomogeneous systems. Moreover, we can apply the 

SMigdal-Kadanoff renormalization[10, 11] scheme to obtain exact solutions of the 
I I partition functions of hierarchical spin models. Among a variety of hierarchical 

'"^ ' spin models, the zero-field g-state Potts model defined on diamond-hierarchical 

^ I lattice has been widely studied. The results revealed from the exact solutions 

O ■ show that the corresponding systems have some nonclassical thermodynamic 

features. They include that the second-order phase transition occurs at finite 
temperature with the cusp-behavior appearing in the specific heat at the critical 
Jv>( , point[8, 9, 12], the distribution of the partition function zeros owns a multifrac- 

^ • tal structure on the complex temperature plane, and the free energy near the 

5^ I critical point has a spatially modulated structure[6, 7]. 

We notice that owing to the infinite iterative steps of decorations in lattice 
construction, there exist a well-defined thermodynamic limit for a hierarchical 
spin model. Thence, the degree of the inhomogeneity is completely specified 
for a hierarchical spin model. If we had performed any finite iterative steps of 
decorations in lattice constrution, it would lead to the lack of thermodynamic 
limit for the model. To improve this while remain the model to be exactly 
solvable, we may study Ising systems defined on decorated lattices. 

A decorated lattice is constructed by implementing bond- or cell- decorations 
to a regular lattice in an iterative way. The geometric complexity of a decorated 
lattice then is characterized by the decoration level n. For a decorated lattice 
of the decoration level n, a primary bond or cell defined on the original regular 
lattice is replaced by an n-bond or cell. An n-bond or cell owns a complicated 
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inherent structure formed by n iterative steps of replacing a bond or cell with 
some basic unit of decoration. In the limit of infinite n, an n-bond or cell 
essentially possesses a fractal lattice as the inherent structure. For this type of 
lattices, the degree of inhomogeneity in the local connectivity is proportional to 
the decoration level n and there exists a well-defined thermodynamic limit for 
any value of n. 

Triangular type decorated Ising models are examples. The analyses on these 
models have been carried out by Plechko[13] and Plechko and Sobolev[14], and 
there exists some non-classical results. In particular, the results of [14] about 
the specific heat of the Ising model on cell-decorated triangular lattices have 
revealed some unusual features of ferromagnetic phase transitions. In the curve 
of the specific heat versus the temperature, as n increases, the critical region is 
narrowing and a round peak appears above the critical point. While the critical 
point marks the set in for the long-range ordering among different n-cells, the 
appearence of the round peak corresponds to the occurrence of the short-range 
ordering inside an n-cell. As n approaches to the limit of infinity, an n-cell has 
the Sierpiiiski gasket as its inherent structure. Then, in approaching the limit, n- 
cells are factorized and become independent in the Ising system. Moreover, the 
critical point approaches the zero temperature, and the round peak temperature 
remains finite. 

Bond-decorated and cell-decorated lattices are different ways of modelling 

inhomogeneous systems. As the decoration level n increases, the order of ram- 
ification remains a constant for a cell-decorated triangular type lattice but it 
increases unboundedly for a bond-decorated square lattice. Furthermore, it was 
pointed out by Gefen et. al. that there is no phase transition at finite temper- 
ature for a fractal lattice with the finite order of ramification[15]. Hence, the 
non-universal features of phase transitions in cell-decorated Ising models may 
be quite different from that of bond-decorations. To obtain a more complete 
picture about the possible features of inhomogeneous systems, we study the 
diamond-type bond-decorated Ising model. Note that the diamond-type bond- 
decoration was also adopted in a widely studied model, the diamond-hierarchical 
g-state Potts model[l — 7]. In this work, we are interested in the following ques- 
tions: {i) How does the degree of the inhomogeneity in a system, characterized 
by the decoration level n, affect the nature of the phase transition? (ii) How do 
the characteristic features of ferromagnetic phase transitions change with the 
lattice structure varying from cell-decorated triangular type lattices to bond- 
decorated square lattices? (Hi) How do the non-universal features obtained 
from the diamond-hierarchical Ising model[l — 9] appear? In particular, can we 
understand better about the cusp-behavior appearing in the specific heat at the 
critical point instead of the well-known logarithmic divergence? 

The analytic expression of the partition function of square Ising model with 
diamond-type bond-decorations can be obtained by referring to the exact result 
of two-dimensional Ising model. The exact solution of the Ising model on a 
square lattice was first given by Onsager[16]. Since then, other alternatives in 
obtaining the exact solution have been developed. Besides of the combinatorial 
method, which was first developed by Kac and Ward[17] and then rigorously 



2 



reformulated by Green and Hurst[18], more recently Plecliko used a nonstan- 
dard method of Grassman multipliers to obtain the analytic expressions of the 
partition functions for a variety of cases[13, 14, 19 — 21]. Based on the frame 
provided by these results, in Section 2 we use the bond-rcnormalization scheme 
to construct the analytic solution of the free energy. In Section 3, we determine 
the critical point of an n-lattice and study how the critical temperature vary 
with n. The condition of the critical point is given by the occurrence of the non- 
analyticity in the free energy. In Section 4, we use the solution of the free energy 
to calculate the internal energy and specific heat in an analytic way. Since the 
free energy can be expressed as the sum of those contributed by the local in- 
teractions among the sub-bonds of an n-bond and the long-range interactions 
among different n-bonds, we also use the same separation in the internal energy 
and specific heat. Thence, in general, we can obtain some insights about the 
nature of the phase transition of an inhomogeneous system by comparing the 
contributions from the local and the long-range interactions. Our results indi- 
cate that the cross over from a finite-decorated system to an infinite-decorated 
system may not be a smooth continuation, we then calculate the critical val- 
ues of the internal energy and specific heat of an infinite-decorated system in 
Section 5. Finally, Section 6 is preserved for the summary and discussion. 

2 free energy 

We study the Ising model defined on an arbitrary n-latticc. A two-dimensional 
square lattice is used as backbones for the construction of an n-lattice. The 
simple square lattice is referred as 0-lattice, and its connecting bonds between 
lattice sites arc named as 0-bonds. For a 0-latticc, the total site-number is 
denoted as rig and the total bond- number as ri;, with rn, = 2ns- Starting from a 
0-bond, we obtain the corresponding n-bond as the result of replacing a bond 
with a diamond iterativcly n times. This n-bond then is specified by the two 
sites connected by the 0-bond. For an n-bond, we have the total site-number 
s(") = 2(4" -I- 2)/3 and bond number 6^"' = 4". The construction of an n-bond 
is schematized in Fig. 1. From the primary structure of a 0-lattice, we replace 
each 0-bond with an n-bond to form an n-lattice. For an n-lattice, the average 
site- and bond-numbers per unit square are A^i"' = 2s("' — 3 and N^"^ = 26^"' 
respectively. Examples of n-latticcs are shown in Fig. 2. 

The general form of the partition function for the Ising model defined on an 
n-lattice reads 



where the sum is over all bond-connected pairs of the n-lattice, and 

the Ising spin takes two possible values ai = ±1. Here we consider uniform 
ferromagnetic couplings characterized by the coupling strength J, and the di- 
mensionless coupling parameter r] is defined as rj = J/ksT. 
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To implement the bond-renormalization scheme in the evaluation of the par- 
tition function, we rewrite Eq. (1) as 



2-(n) 



E 



n 



B 



(") 



(2) 



where the Boltzmann n-factor, B^^^^^^ , is the contribution to the partition func- 
tion from the Boltzmann factors associated with the n-bond specified by the 
two sites, /i and v, the sum is over the spins defined only on the lattice sites 
of the corresponding 0-lattice, and the front factor is added to compensate the 
normalization factor added in Note that the number of the decorated 

sites of an n-bond is s^"-* — 2, and the Ising spins on the decorated sites are 
referred as the inner-spins which couple only to the spins on the same n-bond. 
Thence, the B^i^^^^ of Eq. (2) is defined as the result of taking the sum over the 
inner-spins for the product of all Boltzmann factors associated with the n-bond. 



{<Tfc,fe=l,...,s(")-2} 

(3) 

where the sum is over the - 2) inner-spins, and the front factor is added 
for the normalization of the sum. 

To obtain the explicit form of b'^i^^s^ for an arbitrary n, first we express the 
Boltzmann 0-factor, which takes the form of exp (r?fT^cr„), as 

-^<°!«> = (coshr/) + a^,<Jv (sinhr/) . 
According to the definition of Eq. (3), the Boltzmann 1-factor is given as 



(4) 



(5) 



Using Eq. (4) for exp (ryo-^cr^), after some algebraic manipulations we obtain 
the result of Eq. (5) as 



B 



(1) 



exp(ry^^^) cosh {jl^^^^ + <^ti<^v sinh {jl^^^^ 



with 



exp ('y^^^) = cosh (277) . 
By continuing such construction, we can express Eq. (3) as 

^£1) = ^^"^ (^) [cosh (r,(")) + a^a, sinh (r;^")) 



(6) 
(7) 

(8) 
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for n > 1, where the function i?^"^ (77) is given by 

^^"U^) = n [exp(7?('=))]' \ (9) 

k=l 

and 7]^''^ is the image of r]^''~^^ under the renormalization map defined as 

exp (t?^'^)) = cosh (2ry('=-i)) (10) 

with the initial condition, ry^"^ = 77 for 1 < < n. 

For the Ising system defined on a single n-bond, the partition function of 
Eq. (2) has nt, = 1 and it reduces to 



^(n) ^ 2^(^) f 1 



{fj,,v} 



(11) 



After the substituition of Eq. (8) for , this yields 

Z^") =2^*"'i?(") (r/)cosh . 
Then, the corresponding free energy per bond per fcsT, is given as 

f (") _ 1 i„ vC") 



(12) 



(13) 



where fe*^") is the total bond number contained in an n-bond. This, up to a 
constant, yields 



4n 



(14) 



fe=l 



Note that /^^"^ in the Hmit of infinite n is exactly the free energy density of the 
diamond-hierarchical Ising model. 

Using the result of Eq. (12) for zj^"'\ we can express Eq. (2) in terms of 
Zl"^ as 

(15) 



where Q^^\ referred as the reduced partition function, takes the form of 
Qin) ^ (1)„. ^ J [] [1 + a^a„tanhr?(")] i . 



(16) 



The reduced partition function of an n-latticc given by Eq. (16) has exactly 
the same form as that of simple square lattice when the variable ry^") is changed 
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to rj. Thence, the corresponding free energy can be easily written down accord- 
ing to the formal expression of exact solution of simple square lattice provided 
by Refs. [13, 16, 19,21,22]. Up to a constant, we can express the free energy 
per bond per fc^T as 

= (17) 



where /^^"^ is given by Eq. (14), and fl^^ is the contribution from Q^") and it 



IS given as 



hb 



2nFJo ^Jo 



d0 , 
— In 
27r 



A^"^ (cos 6» + cost/.) 



(18) 



with j4q"^ and A^^' defined as 



(n) 



A- 



(«) _ 



cosh^ (2r?(")) 
cosh^ (r?(")) ' 

2sinh(r/(")) 
cosh^ (r/(")) ' 



(19) 



(20) 



Since f^"'^ is the free energy density of independent n-bonds, /j"'' can be iden- 
tified as the contribution from the interactions among different n-bonds. Thus, 
we refer /^"'■' as the contribution from the local interactions and /j^-* as that 
from the long-range interactions. 

To simplify the expression of fj^^^ , we change the variables, (Ji = {6 + (f)) /2 
and W2 = {0 — (j)) /2, to rewrite Eq. (18) as 

I--K/2 



An) 



coswi costi;2 



with 



o4(") 



A 



(n) 



Then, the integral over uji can be easily performed and the result is 
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{n) 



r^'^ln [1 fl + Vl-K(")2sin2u;') 
Jo [2 ^ ^. 



(21) 
(22) 

(23) 



3 critical point 

The free energy density, /^^"•' , is not completely analytic for the physical region of 

temperature. The non-analyticity arises from the second derivative of f'^^^ with 
respect to temperature at the point k^"^ = 1 for the physical region < k^") < 1. 
This implies that for the Ising model on an n-lattice, the long-range interactions 
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among difFcrcnt n-bonds arc responsible for the occurrence of the Ising phase 
transition at the finite temperature, and the critical point is determined by the 
condition = 1. Thus, we can obtain the bulk critical temperature of the 
ferromagnetic phase transition from the solution of the condition, 

^(") _ 2^(") ^ 0. (24) 

Here, for convenience, we use the notation, =, to denote the equivalence estab- 
lished only at the critical temperature. With the definitions of Aq"' and 
given by Eqs. (19) and (20), we can rewrite the critical condition as the familiar 
form, sinh2r7*^") = 1, or as 

exp (r?(")) = h^°^ (25) 

with 

= + (26) 

for an n-lattice. 

As the consequence of the map of Eq. (10), Eq. (25) yields the critical value 
of the exp (77^""^)) variable as 

exp(r?("-i)) (27) 

for n> 1 with 

1/2 



= (/i(0) + ^(/i(0))2_iJ . (28) 



By applying the recursion relation of Eq. (10) continuously, we can establish 
the result, 

exp (tj^"-'^)) =/iW, (29) 
for 1 < k <n, where h^''^ is determined by the recursion relation 

= (^/i('=-i)+'^(Mfc-i))'-iy^'. (30) 

By taking k = n in Eq. (29) we obtain the reduced critical temperature of 
the Ising transition as 



ksTc 



In I 



' (31) 



for an arbitrary n-lattice. For n = 0, this gives the well-known result, ksTc/J = 
2.269185. As n approaches 00, the /i^") value can be determined by the fixed 
points of the map of Eq. (30). In fact, the value is the repellor of the 

map of Eq. (30), and this repellor locates at h^°°^ = 1.839287. The correspond- 
ing critical temperature is {kBTc/J)^^^ = 1.641018. However, the nature of 
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the phase transition at the rcpellor remains to be answered, and this wiU be 
discussed in Section 5 where the case of infinite decoration-level is discussed. 

The numerical values of {ksTc/J)^ for some different n values are given in 
Table 1 , and the result as a function of n is shown in Fig. 3 with the solid curve 
given by the fitting function {ksTc/J)^ = 1.6410 + (0.6281) exp [- (0.5857) n]. 
This result indicates that the critical temperature decreases exponentially to- 
ward the limiting value 1.6410 as n increases. On the other hand, it has been 
shown [14] that for the cell-decorated triangular Ising model the critical tempera- 
ture as a function of n is given by {ksTc/J)^ = 4/ In {An — 4 In [(n -|- 1) /2] -|- 1} 
for n > 1, and this leads to (kBTc/ J)^ = 4/ln(4n) as n ^ oo. Thus, the de- 
creasing rate of the critical temperature toward the limit of infinite decoration 
level is more faster for the bond-decorated Ising model in comparing with the 
cell-decorated Ising model. 



4 specific heat 

In this section, wc follow the standard procedure, for example as given by Rcf. 
[22], to calculate the specific heat for an n-lattice. The analytic result is ex- 
pressed in terms of the complete elliptic integrals of the first and second kinds. 
The contributions to the specific heat contain those from the local and the 
long-range interactions, and the comparison between these two is given quan- 
titatively. The critical behaviors of the specific heats of different n-lattices are 
also discussed. In particular, we compare the results with the cell-decorated 
case. 

Using the free energy density /^"^ of an n-lattice given by Eq. (17), we 
define the dimensionless internal energy per bond as 



dr] 



(32) 



Due to the separable contributions to the free energy density, we can write the 
internal energy density as the sum of two parts. 
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An) _^ (n) 



+ e 



bb ' 



(33) 



with = a/("V5?? and e<i^^ = df^'^'^ /d^. By using /("^ and of Eqs. (14) 
and (23) , we can perform the direct differentiation to obtain the corresponding 
internal energy density as 
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4*^ y V dr] 



(34) 



and 
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dr] 
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where the cxphcit forms of the functions, dlnA^j^''/dr] and dhiK^'^^ /drj, can be 
obtained by direct differentiations of Eqs. (20) and (22), 



— In^^"^ - /^ CQsh (??(")) 3 sinh (??(")) 
Or] ^ 



and 



with 



d_ 



InK^") = 



I sinh(?^(")) cosh (r?(")) 
'2cosh(2r/(")) 4sinh(2r7(") 



sinh(277(")) cosh (2?7(") 



df] 



2'^ J II [tanh (2r,(^-i)) 




(36) 



(37) 



(38) 



for 1 < A; < n obtained from the map of Eq. (10), and if («(")) is the complete 
eUiptic integral of the first kind defined as 



K 



(«"■') 



Tv/2 



(1-K(")2sin2a;)'^'^ 



(39) 



We notice that K («;(")) is singular at the critical point k^"^ = 1. As k^"^ 
1, it can be shown[22] that 



K 



) = -1" 


T-Tc 


/ T^Tc 


Tc 



(40) 



This logarithmic divergence appears in the specific heat, and it becomes the 
signature of the Ising phase transition in two-dimensions. But, by substituting 
the critical condition, sinh277("' = 1, into Eq. (37), we obtain 



d_ 
df] 



InK^") =0. 



(41) 



As a consequence of Eq. (41), we see that e[^^ of Eq. (35) is a continuous 
function of temperature T even at the critical point, T = Tg. 

To have a qualitative understanding about the varitions of ej;""* and ej;"'' with 

respect to the decoration-level n, we compute the critical values of e^"' and e'^^^ 
for different n values. By using the critical condition of Eq. (29) for an n-lattice, 
we obtain the critical values of Eqs. (36) and (38) as 



(4- 2V2 



and 



dr] 



4 2'= 



(/i(. 



(42) 



(43) 
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for 1 < fc < n. Here and hereafter, we use the subscript c to denote that the 
quantity in the parenthesis is evaluated at the critical point of an n-lattice. 
Then, using Eqs. (41) and (42) we can express the critical value of e["' and e[^' 
given by Eqs. (34) and (35) as 
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fe=i 



dr] 



(44) 



and 



An) 
^bb 



(45) 



with {drj'^^) / drj) ^ given by Eq. (43). 

The numerical results of Eqs. (44) and (45) for some decoration levels n 
are given in Table 1, and the corresponding curves are shown in Fig. 4. These 
results indicate that the local correlations among the 6^"^ bonds of a given n- 
bond are enhenced and the long-range correlations among different n-bonds are 
reduced as n increases. For n > 6, the long-range correlations are almost absent, 
and the system is very closed to that composing of independent n-bonds. 

In calculating the specific heat Cy''' of an n-lattice, we use the dimensionless 
quantity, c^"^ = Cy"^ /ks, defined as 



An) 



V 



dr] 



(46) 



Similar to the internal energy density, we can express c*-"-* as the sum of two 
parts. 



in) _ An) 



C^-' = C, 



(n) 
bb ' 



(47) 



with c, 



in) 



rfdej^^/drj and 



in) 



-rj''de'^^ /drj. 



For the direct differentiation of e["^ of Eq. (34) with respect to rj, we obtain 



dri^'')/dr] 
cosh (ry(")) 



+ 



(tanh (r/"^)) 



» 1 / , , _ /^yr^x I " r / 1 \ / 92^W 

(48) 

where the function, dr]'^^) /drj, is given by Eq. (38), and its further derivative 
yields 



drf 



dn^k) 

dr] 



E 



2'+i 



' i-l 



sinh (477(»-i)) 



W tanh (2ry(-'-i' 



(49) 



for A; > 1. Note that in Eq. (49) the product for the index j in the last 
parenthesis exists only for i > 2. 
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On the other hand, the direct differentiation of e^^^ of Eq. (35) with respect 
to 1] yields 



in) 



7^2 2N^"'^ I \dri'^ 

b 



where the functions, Pq (k^"^) and Qo (k'"^), are defined as 

'-^ («'"') = (^^»'"')-(l;'-'"')^ 



and 



dr] 



ln«;(") 



(50) 
(51) 

(52) 



and E (k^")) is referred as the complete elliptic integral of the second kind and 
defined as 

/ , , X 1/2 



da;(l-K(")2sin2u;) 



(53) 



Note that in obtaining Eq. (50) we have used the well-known formula[22]. 



dK («;(")) 



1 



^1 -(«(")) ^ 



(54) 



We also notice that by further differentiating the results of Eqs. (36) and (37) 
with respect to rj we obtain 



drf 



dr] 



and 



with 



32^(n) 

drf' 
drf 



(n) 



and 



Q2 
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~ sinh^ (ry^")) 


cosh^ (r?(")) ' 


cosh (rj'^^^^ 


3 sinh (??(")) 


sinh (?7(")) 


cosh(?7(")) ' 
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sinh^ (2r/(")) 


cosh^ (2r?(")) 


2cosh(2?7(")) 


4sinh(2??(")) 



sinh(2r7(")) cosh(27?(")) 



(55) 

(56) 
(57) 

(58) 
(59) 

(60) 
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The valucss of c^"-* and c^'^' of Eqs. (48) and (50) can be calculated by using 
the explicit forms of Eqs. (36-38), (49), (55), and (56). The numerical results as 
a function of reduced temperature ksT/J for some n values are shown in Fig. 
5. The peak temperature of c^^\ denoted by {kBT^I J)n, marks the occurrence 
of the local ordering in an n-bond, and the numerical values of {kBTm/ J)^ for 
some decoration levels n are given in Table 1. The plot of {kBTm/ J)„ versus n 
is also shown in the lower part of Fig. 3 with the solid curve given by the fitting 
function {kBTm/J)n = 1.6410- (0.8063) exp [- (0.7144) n]. Thus, because of the 
enhancement of the local correlation for the increase of n, {kBTm/ J)„ increases 

as n increases. Moreover, the peak temperature of c^"' is the critical temperature 
{ksTc/ J)^ which signifies the occurrence of the long-range ordering among the 
Ising spins defined on the primary sites, and the behavior of {kBTd J)^ versus n 
is already given in the upper part of Fig. 3. The results of Fig. 3 indicate that 
the local ordering occurs at the temperature below the critical temperature, and 
this is exactly opposite to the case of cell-decorated triangular Ising model[14]. 

The singular behavior of c^"^ arises from of Eq. (50) in the form of the 
logarithmic divergence given by Eq. (40). Thus, we can express the singular 
behavior of c^"^ as 

(61) 



T^Tc 



^sin g ill 



where A^^^^, referred as the critical amplitude, is given as 



^4^") =- 

sing 



Pf 



(62) 



with Pq («(")) given by Eq. (51). By using Eq. (41), we can rewrite the 
expression of A^^^^ as 



An) _ (d- 



(63) 



By taking the critical value of Eq. (56) with [Q2 (??'"^)] ^ = and [P2 (r?^"))] ^ = 
8, we can use the result of Eq. (43) to obtain 



q2 An+2 

ar?2^ " 2 



n 

, fe=i 



(feW)'-l 



(64) 



Thus, using the reduced critical temperature given by Eq. (31) we can express 
Eq. (63) as 



sing 



In 



9 I " 

"'"'] n 



(/iW)' 



(65) 



The numerical values of Eq. (65) for some n values are given in Table 1. 
The plot of ^sCig versus n is also given in Fig. 6 where the solid curve is given 
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by the fitting function A^^^^^ = (0.2473) cxp [- (0.3018) n]. Since the critical 
amphtude associated with the logrithmic singularity characterizes the width of 
the critical region, our results indicate that as increasing n the critical region is 
narrowing down to the vanish limit. This characteristic feature in the specific 
heat can also be seen in the model on anisotropic regular lattices[13, 14, 20, 22]. 
We also notice that for the cell-decorated triangular Ising model the critical 
amplitude behaves as A^^^ ~ 3~" (nlnn)^ for very large n [14]. Thus, as the 
decoration level n increases the critical amplitude approaches zero less quickly 
for the bond-decorated model in comparing with the cell-decorated model. 



5 infinite decoration-level 

In this section, we discuss the critical properties of the Ising model on an n- 
latticc in the limit of infinite n. First we demonstrate the absence of Ising phase 
transition. This feature is caused mainly by the fragmentation of n-bonds. 
Then, we calculate the critical values of the internal energy and specific heat. 
The specific heat curve near the critical point becomes a cusp, and the exact 
critical value is given. 

The map of Eq. (10) has three fixed points, 77^, given by the solutions of 

rjf =\n [cosh (2r]f)] . (66) 

Among the three fixed points, two are attractors locating at 77°-^ = and 77°^ = 
00 respectively, and one is a repellor locating at ry^ = 0.609378. For the limit 
of infinite decoration-level, the repellor is the critical point of the bulk system. 
Since the critical point is a repellor. the critical behavior of the system may 
be different from that obtained by extrapolating the results of finite-decorated 
systems to the limit of infinite n. To show this, first we give the plots of 77^"' 
versus r/ for 7?, = 0, 2, 4, and 6 in Fig. 7. The results of Fig. 7 indicate that for 
the limit of infinite n the possible 77^^"' values are restricted to the fixed points 
of Eq. (66): 7?(") = rjf for r? < r?^, 7?^") = r?^ for r] = r)}, and 7/(") = rjf for 
Tj > r]^f. This leads to the result that the value of k^"' defined by Eq. (22) has 
only two possible values, 

= (67) 

for 77 ^ T/j, and 



exp(-277p] - [exp (-4?7^)] (68) 



for T] = rfp The latter gives the numerical value k*^""' = 0.912622. Thence, 
the logarithmic divergence of the Ising transition, caused by the divergence of 
the elliptic integral K (k'^"^) of Eq. (39) at the point k*^"^ = 1, is absent for 
the infinite n limit. This indicates that the cross over from a finite-decorated 
system to an infinite-decorated system may not be a smooth continuation. 
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To exhibit the critical properties of the infinite-decorated system more ex- 
phcitly, we then calculate the critical values of the internal energy and specific 
heat. Since the critical point of the infinite-decorated system is the repellor 
of the map of Eq. (10), we have the following identities for the derivatives 
evaluated at the critical point j] = rfj: 



drj' 



(fe) 



and 



(69) 



(70) 



Here, we use the notation, = , to denote the equivalence established at the 
repellor. On the other hand, by successively applying the chain rule we have 



dr] J--^ ■ 

Then, as a result of Eq. (69) we obtain 



(71) 



(72) 



dr] 

with = {dr]^^^/dr])^^^, = 1.678574. 

For the second derivative, d'^rj^^'' / drf , we first perform the direct differenti- 
ation of Eq. (71) and then use the chain rule to obtain 



with 



drj^ 



dr] 



32^(2) 

dr](- 
dr](-^ 



,(2)\ n / /)2 (fe) \ 



(73) 

(74) 



for a < k. Then, by using the identities of Eqs. (69) and (70) we can express 
Eq. (73) as 



dr]"^ 



2(fe-l)' 



which yields 



dvp 



= d(2) 



fe-i 



1 



- 1 



(75) 



(76) 



with ^ (5277(1) = 1.182391. 
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To calculate the internal energy density and specific heat at the repellor, we 
take ?7 = 77^., Eq. (72) for {drj^'^y drj) , and Eq. (76) for [d^f]^''^/ drf) . 

Then, for the internal energy density of Eqs. (34) and (35) we obtain the 
numerical results as 



^ - (0.543689) (a(^))" - a^^) (^y^^) , (77) 



and 

An) 



9' 



= - (0.414742) (a(i))", (78) 

with a(i) = d(i)/4 < 1. Similarly, for the specific heat of Eqs. (48) and (50) at 
the repellor we obtain the numerical results as 

4"^ = (0.261574) (a(2))"+(0.201894) 5(1) 1^- -^-L_^ (a'^)) "+(0.371342) 5(1) 

(79) 

and 

^ - (0.115082) (a(2))" + (0.154011)^(1) - j (a*^))" , (80) 

with a(2) = '74 < 1, 5^') = rf('V {{d^^^f - d^^)),9^^^ = a(2) (l - (a(2))") / (l 
and g^^) = a(i) (l - (a(i))") / (l - a(i)). 

For the limit of infinite decoration-level, we have (a«)" ^ 0, (a(2))" ^ 0, 
and 1/ ((id))" ^ as n — > 00. Then, in this limit Eqs. (77) and (79) reduce 
to 

^^"^ = - (1^) = -0-723079, (81) 



and 

a(2) \ / aW 



4°°) = (0.371342)^(1) 



(2)/ Vl-a(i 



0.639852, (82) 



while the long-range parts vanish, e|;^) = and c[^) = 

From these results we may conclude that the system of an n-lattice reduces 
to that composing of independent n-bonds as n — > 00. In this limit, our results 
also show that the critical specific heat becomes a cusp with the height c(°°) = 
0.639852 at the critical point. Thus, the corresponding critical exponent of the 
specific heat, a, is the same as that of the diamond- hierarchical Ising model, 
a= -0.67652 [9, 23]. 
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6 summary and discussion 



We use the Ising model defined on square lattices with diamond-type bond- 
decorations to study the ferromagnetic phase transition of an inhomogeneous 
system. In this model, the decoration level n associated with a decorated lattice 
can be used as a parameter to characterize the degree of the inhomogeneity of 
a system. Based on the cixact solutions of the two-dimensional Ising mobel, 
we employ the bond-renormalization scheme to obtain the analytic form of the 
free energy. The free energy and its derived quantities, the internal energy 
and specific heat, are expressed as the sum of those contributed by the local 
interactions among the sub-bonds of an n-bond and the long-range interactions 
among different n-bonds. Then, the change in the nature of the phase transition 
as n increases can be viewed as the reflection of the competetion between the 
local and the long-range interactions. 

In general, as n increases, the long-range interaction is reduced and the local 
interaction is enhenced. The logarithmic singularity of the specific heat at the 
critical point is caused by the long-range interaction. Along with the reduction 
of the long-range interaction as n gets larger, we have the critical temperature 
lowering down and the width of the critical region becoming narrower. By fitting 
the data numerically, we obtain the critical temperature changes with the deco- 
ration level n as (fcsTc/ J)„ = 1.6410+(0.6281) exp [- (0.5857) n] and the critical 
amplitude of the specific heat changes as A^^^^g = (0.2465) 4" exp [— (1.6530) n]. 
Moreover, because of the enhancement of the local interactions, the occur- 
ring temperature of the local ordering increases with the form, {kBTm/ J)n = 
1.6410-(0.1081) exp [- (0.4098) n], as n increases. Both {kBTc/J)^ and (kBT,n/J) 
approach to the same limit 1.6410 as n approaches the infinite limit. Thus, the 
distance between {UbTc/J)^ and {kBTm/J)^ decreases with increasing n. 

The cross over from a finite-decorated system to an infinite-decorated sys- 
tem is not a smooth continuation. For the infinite decoration-level, the critical 
temperature fcsTc/J is 1.64101(8) which is exactly the location of the repellor 
of the renormalization map of the variable r]^'^\ In this limit, the system of an 
n-lattice reduces to that of independent n-bonds, and the specific heat becomes 
a cusp with the value 0.63985(0) at the critical point instead of the logarithmic 
divergence. 

It is interesting to compare the results we obtain with those from the cell- 
decorated triangular Ising model given in Ref. [14]. The critical region decreases 
with increasing n for both cases. But, there exists some different features as 
the foUowings: (i) The decreasing speed of the critical amplitude of the specific 
heat with increasing n is slower for the bond-decoration case in comparing with 
that of the cell-decoration, (m) For the case of bond-decorations, the occurring 
temperature of the local ordering {kBTm/ J)^ of an n-lattice is smaller than 
the critical temperature (fc^Tc/J)^ and the distance between {kBTm/ J)n 
{kBTc/ J)^ decreases with increasing n. However, for the case of cell-decorations, 
{kBTm/ J)n is larger than {kBTc/ J)^ and {kBTm/ J)^ is freezed for n > 5. {iii) 
There is a finite peak arising to the right of the critical point for the specific 
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heat due to the local ordering in the cell-decoration case, but such a peak is not 
observed in the bond-decoration case. 

We may use the idea of the order of ramification to explain the differences 
described in the above. To characterize the geometric difference between bond- 
decorated and cell-decorated lattices, we introduce the quantities, i?iand i?2, 
with i?i for the cutting-bond number of isolating an n-bond or n-cell and R2 
for the average of the total bond number contained in an n-bond or n-cell. 
We notice that the Ri value signifies the ability of establishing the long-range 
correlations among different n-bonds or n-cells and the interactions among R2 
bonds are referred as the local interactions. 

For the cell-decorated triangular Ising model with decoration-level n, we 
have Ri = 12 and R2 = 3"+^. Because of Ri <^ R2 for n > 3, the long-range 
correlation is very hard to establish, and this leads to the higher temperature for 
the occurrence of the local ordering in comparing with that of the long-range 
ordering. For sufficiently large n with R1/R2 — 0, an n-cell almost become 
a closed subsystem, and thence a round peak, contributed by the subsystem, 
appears above the critical point in the curve of the specific heat. 

For the diamond-type bond-decorated Ising model discussed in this work, 
we have i?i =6-2" and R2 =4". The values of Ri and R2 are comparable 
for n < 5, and the long-range effect is almost absent for n > 6. Thence, the 
round peak in the specific heat, which appears in the case of cell-decorations, 
is absent. Moreover, the coordination numbers of the decorated sites in an n- 
bond are highly inhomogeneous, and this is mainly responsible for the lower 
occurring temperature of the local ordering in comparing with that of the long- 
range ordering. 

Thus, the non-universal features of phase transition for the diamond-type 
bond-decorated Ising model is quite different from that for the cell-decorated 
triangular Ising model. 
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Table 1. The numerical results of the crical temperature (k-BTc/ J)^, the 
occurring temperature of the local ordering (kBTm/J)^, the critical values of 
the internal energy density (e^"^)^, the internal energy density due to the local 

interactions (^ej,"^^ ) and the internal energy density due to the long-range in- 
teractions (^<^ifa^ j ) and the critical amplitude at the logarithmic divergence in 
the specific heat A^^^ . 
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Figure 1: Examples of n-bonds specified by the two lattice sites /i and v for 
(a) n = 0, (b) n = 1, and (c) n = 2. Note that the Ising spins except and 
Uv are referred as the inner-spins. 



Figure 2: Examples of n-lattices for (a) n = 0, (b) n = 1, and (c) n = 2. 
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Figure 3: The critical temperature (ksTc/ J)^^ and the peak temperature of the 
specific heat {kBTm/ J)n caused by the local interactios of an n-lattice denoted 
by black dots and black triangles respectively. The continuous curves are given 
by the fitting functions as (/.-sTc/J),^ = 1.6410 + (0.6281) exp [- (0.5857) n] and 
{kBTm/J)n = 1.6410 - (0.8063) exp [- (0.7144) n] respectively. 
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Figure 4: The critical values of e*^"^ (triangles), (black dots), and ej 
(squares) as a function of decoration-level n. 
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Figure 5: The numerieal results of cj^ (with round peaks) and cj^ (with sharp 
peaks) as a function of reduced temperature ksT/J for n = 2 , 4, and 6. 
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Figure 6: The critical amplitude A^^ng of the specific heat as a function of 

decoration-level n with the solid curve given by the fitting function A^^^g = 
(0.2473) exp [- (0.3018) n]. 
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Figure 7: The variable as a function of for n = 0, 2, 4, and 6. Note 
that the vertical line corresponds to the location of the repellor of the map 
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